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1 | INTRODUCTION 


W.-Y. Gu® | S.-X.Luo* | Q.-K. Shi? | L.-S.Hou? | 


Summary 

Porcine deltacoronavirus (PDCoV) is a recently identified coronavirus in the genus 
Deltacoronavirus that can cause enteric disease with clinical signs including diar- 
rhoea, vomiting, dehydration and mortality in neonatal piglets. Although evidence of 
the prevalence of PDCoV in China is accumulating, little published information 
about Chinese PDCoV isolates is available. In this study, we investigated the pres- 
ence of PDCoV in 49 faecal/intestinal samples from piglets with diarrhoea on differ- 
ent farms in Hebei province. Five samples (10.2%) were positive for PDCoV, but no 
coinfection of PDCoV with other enteropathogens was observed. A PDCoV strain 
named HB-BD was successfully isolated from the intestinal contents of a diarrhoeic 
piglet and serially propagated in swine testicular (ST) cells for >40 passages. The 
complete genome of the HB-BD strain was sequenced and analysed. Genomic anal- 
ysis showed that the HB-BD strain had a closer relationship with Chinese strains 
than those from other countries and was grouped within the Chinese PDCoV clus- 
ter. The results of this study will be valuable for further research of PDCoV genetic 
evolution and development of effective diagnostic reagents, assays and potential 


vaccines against newly emerged PDCoV strains. 
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with lower mortality in piglets (30%-40%) than typical PEDV infec- 


Porcine deltacoronavirus (PDCoV) is a novel member of the genus 
Deltacoronavirus of the family Coronaviridae. The disease caused by 
this virus is characterized by severe watery diarrhoea and vomiting 
in neonatal piglets. It is symptomatically indistinguishable from diar- 
rhoea caused by porcine epidemic diarrhoea virus (PEDV) and trans- 
missible gastroenteritis virus (TGEV) (Li et al., 2014; Song et al., 
2015; Wang, Byrum, & Zhang, 2014), but presents more mildly and 
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tion (Hu, Jung, et al., 2015). Since its initial isolation in the United 
States in February 2014 (Li et al., 2014; McCluskey et al., 2016; 
Marthaler et al., 2014), PDCoV has been identified in Canada, South 
Korea, China and Thailand. Molecular surveillance of diarrhoeal sam- 
ples of pigs from the United States indicated a 30% infection of 
PDCoV, which is considered a common viral pathogen of pigs 
(Marthaler et al., 2014; Zhang, 2016; Jung, Hu, & Saif, 2016). How- 
ever, owing to the fact that the PDCoV is a recently emerging viral 
pathogen, there is limited knowledge about the virus (Thachil, Ger- 


ber, Xiao, Huang, & Opriessnig, 2015). There are presently no 
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effective treatments or vaccines available to control PDCoV (Hu, 
Jung, et al., 2015). 


The PDCoV is an enveloped, positive-sense, single-stranded RNA 
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virus. The full genome of PDCoV is approximately 25.4 kb in length 
(Woo etal., 2012), encodes two polymerase proteins (ORF1a/b) 
gene, four major structural proteins for the spike (S), envelope (E), 
membrane (M) and nucleocapsid (N) proteins, as well as non-struc- 
tural protein 6 (NS6) and non-structural protein 7 (NS7) (Lee & Lee, 
2014; Li et al., 2014). The functions of PDCoV individual proteins 
have not been elucidated. However, among the structural proteins 
of the other coronaviruses, the S protein is known to play pivotal 
roles in interacting with the cellular receptor to mediate viral entry 
and is an antigenic target for neutralizing antibodies (Park, Song, Ha, 
& Park, 2007; Sato et al., 2011). The M protein is a transmembrane 
protein that is responsible for the transport of nutrients across the 
membrane. In addition, the M protein plays an important role in viral 
assembly, as well as the induction of virus neutralization antibodies 
(Fan, Zuo, Li, & Pei, 2012). The coronavirus N protein binds viral 
RNA and highly conserved. It is involved in several of the biological 
activities of the virus (Chang et al., 2006; Molenkamp & Spaan, 
1997). The S, M and N protein genes have been targeted for the 
development of virological and serological diagnostic assays for 
PDCoV (Chen, Gauger, et al., 2015; Ma et al., 2015; Marthaler et al., 
2014; Song et al., 2015; Su et al., 2016; Wang et al., 2014). There- 
fore, it is important to understand the molecular characteristics of 
the prevalent strains and its relevance to the development of diag- 
nostic reagents, assays and potential vaccines against emergent 
PDCoV strains. Although some publications have reported the detec- 
tion of PDCoV RNA in domestic pigs in mainland China (Chen, Zhu, 
et al., 2015; Song et al., 2015; Wang, Yue, Fang, & Huang, 2015), a 
few have reported on the isolation of PDCoV or the characteristics 
of those isolates. 

In this study, five of 49 samples collected from different swine 
farms reporting cases of diarrhoea in Hebei province were PDCoV 
positive. A PDCoV strain, HB-BD, was successfully isolated and seri- 
ally propagated in swine testicular (ST) cells. The complete genome 


of the HB-BD strain was sequenced and analysed. 


2 | MATERIALS AND METHODS 


2.1 | Sample collection and testing 


From October 2015 to March 2016, faecal/intestinal samples 
(n = 49) from piglets with diarrhoea on different farms in Hebei pro- 
vince of China were submitted to the Veterinary Preventive Medicine 
Laboratory, College of Veterinary Medicine, Hebei Agricultural 
University, Baoding, China. To detect the PDCoV genome in the col- 
lected samples, viral RNAs were extracted from samples using the 
TRIzol reagent (Invitrogen, Carlsbad, CA, USA) and reverse transcrip- 
tion was performed using the PrimeScript™ 1st strand cDNA Synthe- 
sis Kit (TaKaRa Biotechnology, Dalian, China) with an oligo dT primer. 
RT-PCR primers (F: 5’/-TAACTCCGCCATCAAAC-3’ and R: 5’- 
CCACTTCCACGCTCCT-3’) targeting the N gene were used with the 
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following reaction conditions: 94°C for 5 min, followed by 35 cycles 
at 94°C for 30 s, 50°C for 30 s and 72°C for 30 s, with a final exten- 
sion at 72°C for 10 min. In addition, molecular detection of the other 
three diarrhoea-related enteric viruses, PEDV, TGEV and porcine 
rotavirus (PRoV), was performed using previously reported methods 
for further evaluation of possible coinfection status with PDCoV (Hu, 
Li, et al., 2015; Jeong et al., 2009; Sun, Leng, Zhai, Chen, & Song, 
2014). 


2.2 | Virus isolation, propagation and titration 


Isolation of PDCoV was attempted on samples that were positive 
for PDCoV alone. The isolated strain was used for sequence analysis 
and phylogenetic analysis, the details of which are presented in 
Table 1. To get the inoculum for virus isolation, 1 ml of the original 
intestinal sample was diluted with phosphate-buffered saline 
(0.01M, pH7.4), vortexed and centrifuged at 4,200 g for 10 min at 
4°C. The supernatants were separated and filtered through a 
0.22 um filter and then used as the inoculum for virus isolation in 
cell culture. 

The PDCoV was isolated in the ST cell line (American Type Cul- 
ture Collection No. CRL1746). The growth medium for ST cells was 
advanced Dulbecco’s modified Eagle’s medium (DMEM) supple- 
mented with 10% heat-inactivated foetal bovine serum (Zhejiang 
Tianhang Biotechnology Co. Ltd. Hangzhou, China) and 1% peni- 
cillin-streptomycin solution, and the maintenance medium for 
PDCoV propagation was DMEM supplemented with 1% penicillin— 
streptomycin solution. 

For the first inoculation, cells were cultured in a T25 flask. When 
the ST cell reached 80% confluence in the flask, cells were washed 
twice with maintenance medium and were used for virus inoculation. 
A volume of 900 ul of maintenance medium together with 100 ul of 
the filtered sample was added to the flask. After incubation of the 
virus for 60 min at 37°C in 5% CQO, incubator, the cells were 
washed twice, and 10 ml maintenance medium with 1% pancreatin 
was added to each flask. The cells were then cultured continuously 
at 37°C in 5% CO, incubator and were observed daily for the cyto- 
pathic effect (CPE). When an obvious CPE was observed in around 
90% of the cell monolayers 2-3 days after inoculation, the flasks 
were frozen at —80°C and thawed three times. The supernatants 
and cells were harvested together and stored at —80°C. These sam- 
ples were used as seed stocks for subsequent passage and for the 
detection of PDCoV and other enteric viruses. The isolated PDCoV 
strain, HB-BD, was plaque purified as previously described (Hu, Jung, 
et al., 2015). 

For PDCoV serial propagation, 80% confluent ST cells in a T25 
flask were washed twice with maintenance medium, and then, 1 ml 
of the virus inoculum containing 1% pancreatin was added to 
the flask. The cells were cultured at 37°C in 5% CQ», incubator. 
After 2-3 days, 80% of the virus-infected cells showed CPE, and 
the cultures were collected after three frozen-thawed cycles. The 
supernatants were stored at —80°C and used for PDCoV propaga- 


tion. 
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TABLE 1 Oligonucleotide primers used 
for amplification of the complete genome 
of PDCoV strain HB-BD by RT-PCR 


Primer 
PDCoV1F 
PDCoV1R 
PDCoV2F 
PDCoV2R 
PDCoV3F 
PDCoV3R 
PDCoV4F 
PDCoV4R 
PDCoV5F 
PDCoV5R 
PDCoV6F 
PDCoV6R 
PDCoV7F 
PDCoV7R 
PDCoV8F 
PDCoV8R 
PDCoV9F 
PDCoV9R 
PDCoV10F 
PDCoV10R 
PDCoV11F 
PDCoV11R 
PDCoV12F 
PDCoV12R 
SF1 

SR1 

SF2 

SR2 

SF3 

SR3 

SF4 

SR4 

SF5 

SR5 
PDCoV13F 
PDCoV13R 
MF 

MR 
PDCoV14F 
PDCoV14R 
NF 

NR 
PDCoV15F 
PDCoV15R 


Sequence(5’-3’) 


ACATGGGGACTAAAGATAAAAATTATAG 
GATACTTCATAACTCTGGCAATC 


TGATGTTCTGCTAGCCTT 
CGAGTGTCAGAGGTGTGT 
CACTGATGTAGGCGATGA 
CACGACTTTACGAGGATGA 
TCCATTTGGACCCACCTC 
TAGCCTGCTGACTAAGACG 


AGTCAGCAGGCTATACGTGTGA 
GAATGTTGTCTACTGCCCACGC 


GGAGGCGGTTCACAGTTGTA 


GGTGGAAACCGTAACATTGCTG 


TGGCAGTTAAGATGTCCC 
TGTAGCATTCCTCCTCGA 
CTAACTGCGCTCGGTTTA 
GGTAGAATCGCTGGCTTT 
CTGTGGCAGGAGTGTCTA 
GAAGTTTAATGAAGCGTTG 
AGTCATACGCTACCGCAACC 
GCCATCGGCAACTCCTACT 
TGGCATGATTGTGGTGCAG 
AACAGCTGTGTAGTTGGCAG 
CAACCGCACTAACTTACCTGT 
TT CTGGTGGCCTCACAAGAA 
TATTATCTCGGCTCGTGAG 
AGTGTTATGAGTGTATCGG 
TACTTTACCTATTACGGT 
TGTGATAGCACCGACAACG 
CACAGGTGAGCTTTATGC 
AGAGCCAGTATACATTGCC 
TCTAGAGACATGGCCATCG 
ACTGGTAGAGTATAAGTTG 
TTTTCATGCATGCAGTGCT 
TTAACACAAGCTAAGCAAG 
TGTCATCTGCATTGGTGTGGC 
AGCCTCCTTGGAGGATGACG 
TACTCCAAGCCGAACCCCGT 
GCTGCAAGTGGCAGTTGCA 
CACTTGCAGCTGCGAGATTTA 
GGACTACTGGTGCAGCCATG 
CATCGCTCCAAGTCATTCTT 
CATAGGTTGATGTCTACGCT 


GAGATAAAGCAGGAATCAGCAGC 


GCTCCATCCCCCCTATAAGCC 


PDCoV, porcine deltacoronavirus. 


Nucleotide position? 


1-28 
1,558-1,580 
1,485-1,502 
3,329-3,346 
3,282-3,299 
4,844-4,862 
4,/27-4,744 
6,394-6,412 
6,400-6,421 
8,163-8,184 
7,996-8,015 
9,/81-9,802 
9,710-9,/27 
11,647-11,664 
11,539-11,556 
13,333-13,350 
13,060-13,077 
14,661-14,679 
14,623-14,642 
16,106-16,124 
16,056-16,074 
17,718-17,737 
17,654-17,674 
19,306-19,325 
19241-19259 
20,049—20,067 
19,963-19,980 
20,809—20,827 
20,710—20,727 
21,508—21,526 
21,419-21,437 
22,208—22,226 
22,110—22,128 
22,895-22,913 
22857-22877 
23101-23120 
22,972-22,991 
23,698—23,/16 
23707-23727 
23999-24018 
23,942-23,961 
25,022—25,041 
25002-25024 
25400-25420 


“Nucleotide position is numbered based on the PDCoV/NH strain (KU981059). 


Product size (bp) 
1580 


1862 


1581 


1686 


1785 


1807 


1956 


1812 


1620 


1502 


1682 


1672 


827 


867 


817 


808 


804 


264 


744 


312 


1110 


419 
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Viral titre was measured using 50% tissue culture infectious dose 
(TCIDs9) assays on ST cells in a 96-well plate, as previously 
described (Hu, Jung, et al., 2015). 


2.3. | Immuno fluorescence assay 


PDCoV-infected ST cells were fixed with 4% paraformaldehyde at 
4°C for 30 min then washed three times with PBS and permeabilized 
with 0.2% Triton X-100 for 10 min at room temperature. After 
blocked with 5% bovine serum albumin (BSA) at 37°C for 1 hr, the 
cells were incubated with PDCoV polyclonal antibodies which were 
produced from PDCoV HB-BD-immunized BALB/c mice. After 1 hr 
incubation at 37°C, the cells were washed three times with PBS and 
incubated with FITC-labelled goat anti-mouse secondary antibody 
for 1 hr. The stained cells were examined using a fluorescence 


microscope. 


2.4 | Sequencing and phylogenetic analysis of the 
complete genome of the HB-BD strain 


The 5th passage of PDCoV HB-BD cell culture sample was cen- 
trifuged at 4,000 g for 5 min. A volume of 250 ul of the clarified 
supernatants was used to extract viral RNA using the TRIzol reagent 
(Invitrogen) according to the manufacturer’s instructions. The RNA 
samples were stored at —80°C and used as the template for RT- 
PCR. 

To characterize PDCoV HB-BD genome, a set of overlapping pri- 
mers were designed based on the PDCoV strain, NH (GenBank 
accession no. KU981059), to amplify the nearly complete PDCoV 
genome (Table 1). The primers were synthesized by Beijing Sun- 
biotech Co. Ltd. (Beijing, China). Reverse transcription was performed 
using the PrimeScript™ ‘1st strand cDNA Synthesis Kit (TaKaRa 
Biotechnology, Dalian, China) and carried out using an oligo dT pri- 
mer. The cDNA was amplified by PCR using Taq DNA Polymerase 
(TaKaRa, Dalian, China). The PCR mixture (20 wl) contained 4 ul 
cDNA, 2 ul of 10 x PCR Buffer, 3 ul of 2.5 mM dNTP, 1 ul forward 
primer and reverse primer, respectively, 0.5 ul of Taq polymerase (5 
U/ul) and 8.5 ul of HzO. Fragments of the HB-BD strain were ampli- 
fied under the following conditions: denaturation at 94°C for 5 min, 
35 cycles of 94°C x 30s, 55°C x 305 and 72°C x 2 min, and a 
final extension of 72°C for 10 min. 

The PCR product was cloned into the pMD19-T cloning vector 
(TaKaRa). The positive recombinant plasmids were verified by restric- 
tion enzyme digestion and were submitted to a commercial company 
(Sangon Biotech, Shanghai, China) for sequencing. The raw genomic 
sequence fragments were imported to SeqMan in DNASTAR (DNAS- 
TAR, Inc., Madison, WI, USA) for assembly and annotation. 

Sequence alignment analysis was processed using the Clustal W 
program in the DNAStar software, and the phylogenetic tree was 
constructed by the neighbour-joining method with 1000 bootstrap 
replications using the MEGA 7.0.14 software. The PDCoV/HB-BD 
strain and the reference strains used for the sequence analysis are 


presented in Table 2. 


TABLE 2 The HB-BD isolate strain and other reference strains of 
PDCoV and other coronaviruses used for sequence alignment and 
phylogenetic analysis 


Location and 


temporal Accession 
Strain information no. 
PDCoV/HB-BD-complete genome — China MF948005 
PDCoV/HB-BD-M gene China KY129985 
PDCoV/HB-BD-N gene China KY129986 
PDCoV/HB-BD-S gene China MF037204 
CHN-JS-2014 China KP757892 
CH-HB-2014 China KP757891 
CH/SXD1/2015 China KT021234 
PDCoV/NH China KU981059 
CHN-Tianjin China KY065120 
CH/Sichuan/S27/2012 China KT 266822 
CHN-AH-2004 China KP757890 
HKU15-155 China JQ065043 
HKU15-44 China JQ065042 
PDCoV/CHJXNI2/2015 China KR131621 
OH11846 USA KT381613 
OH1978 USA KJ462462 
PA3148 USA KJ584358 
IL2768 USA KJ584355 
NE3579 USA KJ584359 
M16148 USA KJ620016 
IN2847 USA KJ569769 
SD3424 USA KJ584356 
KY4813 USA KJ584357 
OhioCVM 1/2014 USA KJ769231 
8734/USA-IA/2014 USA KJ567050 
PDCoV/USA/Illinois133/2014 USA KJ601777 
PDCoV/USA/llinois136/2014 USA KJ601779 
PDCoV/USA/llinois134/2014 USA KJ601778 
PDCoV/USA/Illinois121/2014 USA KJ481931 
KNU14-04 Korea KM820765 
PDCoV/Swine/Thailand/S5011/ Thailand KU051641 

2015 
PDCoV/Swine/Thailand/S50O15L/_ —_—s‘ Thailand KU051649 
2015 

PEDV/CH/JX-1/2013 China KF760557 
PEDV/CV777 Belgium AF353511 
PEDV/OH1414 USA KJ408801 
TGEV H16 China FJ755618 
TGEV-virulent-Purdue USA DQ811789 


PDCoV, porcine deltacoronavirus. 


2.5 | Nucleotide sequence accession number 


The nucleotide sequences of the complete M, N, S gene and com- 


plete genome of the PDCoV strain, HB-BD, were deposited in the 
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GenBank database with the accession nos. KY129985, KY129986, 
MFO37204 and MF948005, respectively. 


3 | RESULTS 


3.1 | PDCoV detection 


Clinical samples collected from diarrhoeic piglets from different 
farms in Hebei, China, between October 2015 and March 2016 were 
submitted for pathogen detection and isolation. Of the 49 porcine 
faecal/intestinal samples examined, five (10.2%) were PDCoV posi- 
tive; 31 (63.3%) were PEDV positive; two (4.1%) were PRoV posi- 
tive; and none of the samples were positive for TGEV. Although 
previous studies report coinfection of pigs with PDCoV and PEDV 
(Song et al., 2015; Zhai et al., 2016; Zhang, 2016), none of the sam- 
ples in the present study showed dual infection with PDCoV and 
PEDV. Furthermore, the five PDCoV-positive samples all showed 
single infection with PDCoV. 


3.2 | Viral isolation and propagation in ST cell 
monolayers 


The five PDCoV-positive samples were inoculated on ST cell mono- 
layers, respectively. However, only PDCoV strain HB-BD-inoculated 
cell monolayers showed a visible CPE that consisted of enlarged, 
rounded, densely granular cells that occurred either singly or in clus- 
ters, with evidence of cell shrinkage and detachment (Figure 1). The 
purification of PDCoV HB-BD was performed as_ previously 
described (Hu, Jung, et al., 2015; Hu, Li, et al., 2015). To confirm 
whether PDCoV alone was replicated in ST cells, the cultures were 
evaluated for PDCoV, PEDV, TGEV and PRoV, respectively. Results 
showed that only PDCoV was detected from the HB-BD inoculated 
cell cultures that were all negative for other enteric viruses. This 
demonstrated that the PDCoV HB-BD strain was successfully iso- 
lated from the ST cells. The propagation of PDCoV HB-BD in ST 
cells was also confirmed by IFA staining using mouse anti-PDCoV 
antibodies (Figure 2). To date, the isolate has been serially propa- 


gated in ST cells for >40 passages. 


3.3 | Virus titration performed with the TCIDso 
protocol 


During serial passages, the PDCoV infectious titres were determined 
using a TClIDsq assay on ST cells. Four passages of the isolate (Ps, 
Py, P13 and P19) were detected. The results showed that the virus 
titres for these four passages were 10°° TCIDso/ml, 10°? TCIDso/ 
ml, 10°? TCIDs59/ml and 10°" TCIDs9/ml, respectively. 


3.4 | Sequence alignment and Phylogenetic analysis 
of PDCoV strain HB-BD genome 


The complete genome sequence of PDCoV HB-BD has been 


sequenced and deposited in GenBank under the accession number 


of MF948005. Sequencing results showed that the complete gen- 
ome sequence of PDCoV HB-BD was 25,420 nucleotides (nt) in 
length excluding the 3’ poly (A) tail. The gene order of its genomic 
structure was 5’-UTR-ORF1-S-E-M-NSP6-N-NSP7-3/UTR, and the 
nucleotides numbers of these parts were 540, 18803, 3480, 252, 
654, 285, 1029, 603 and 392, respectively. 

The complete genome sequence of PDCoV HB-BD shared 
97.2%-99.5% nucleotide identity with the other 28 PDCoV refer- 
ence strains (including American and individual Asian strains) avail- 
able in GenBank and had the highest nucleotide identity (99.5%) 
with PDCoV/NH. Both PDCoV HB-BD and PDCoV/NH had a 
unique base insertion in the 5’-UTR (150C151). By comparing S gene 
with those PDCoV reference strains, a 3-nt (AAT, from 154 to 156) 
deletion was observed in the S gene of the HB-BD strain, which was 
also present in other Chinese PDCoV strains, with the exception of 
PDCoV strains, CHN-AH-2004 and HKU15-44. This unique feature 
in the S gene can be used as a genetic marker to discriminate 
PDCoV strains in China from those in the United States, South 
Korea and Thailand (Chen, Zhu, et al., 2015). Further analysis of the 
S gene showed that the HB-BD strain reported in the present study 
shared 95.8%-99.1% nt identity with the other 28 PDCoV reference 
strains. When compared with the Chinese PDCoV strains, the HB- 
BD strain showed 97.9%-99.1% nt identity and 97.2%-99.2% 
deduced amino acid (aa) identity. The HB-BD strain shared the high- 
est nt identity (99.1%) and deduced aa identity (99.2%) with 
PDCoV/NH, which was isolated from Heilongjiang province. Both 
HB-BD and PDCoV/NH isolates had two amino acid mutations 
(L45H and Y123H). Furthermore, HB-BD had four further amino acid 
mutations (H149Y, R888T, A 894T and G910C). 

There were no nucleotide insertions or deletions in the M gene 
ORF of the HB-BD isolate, and the sequence was deposited in Gen- 
Bank under the accession number KY129985. In comparison with 
the foreign strains, we discovered that all Chinese strains, including 
the HB-BD strain, had three nucleotide mutations (G72A, C174T 
and C459T) that were similar to those of the strains from Thailand. 
In addition, the HB-BD strain had a unique nucleotide mutation 
(C63T). However, these nucleotide mutations did not lead to 
changes in the predicted amino acid sequences of the PDCoV iso- 
lates. 

The multiple sequence alignment results of the M gene showed 
that the HB-BD isolate shared high nucleotide homology (99.1%— 
99.4%) with other Chinese strains and had the highest nucleotide 
homology (99.4%) with four Chinese PDCoV strains (HKU15-155, 
JS-2014, SXD1-2015 and NH) and the lowest nucleotide identity 
(98.6%) with two strains from Thailand, $5011 and $5015L. Further- 
more, the HB-BD strain shared 98.9%-99.2% and 99.1% nucleotide 
identities with reference U.S. strains and Korean strains, respectively. 
Moreover, the HB-BD isolate exhibited 100% aa identities with all 
reference sequences except that for the HKU15-155 strain, because 
that strain had one aa substitution (A83P). 

The N gene sequence of all 29 PDCoV strains was identified as 
1029 nt in size that encoded a protein of 342 aa. Sequence align- 


ment results suggested that there was no deletion or insertion in the 
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FIGURE 1. Cytopathic effects of PDCoV isolate in incubated ST cells. ST cells were incubated with the 28th passage of PDCoV HB-BD. (a), 
(b) and (c): mock-inoculated ST cells at p.i. day 1, 2 and 3, respectively. (d), (e) and (f) :HB-BD-inoculated ST cells at p.i. day 1, day 2 and day 3, 


respectively 


FIGURE 2 _ Detection of PDCoV isolate HB-BD (28th passages) in ST cells by IF straining using mouse anti-PDCoV antiserum. (a) IF 
straining of HB-BD-infected ST cells; (b) IF straining of mock-incubated ST cells. ST cells were fixed at p.i. day 1 


N gene regions. The PDCoV isolate, HB-BD, shared 97.8%-99.4% nt 
identity and 98.8%-99.4% aa identity with reference PDCoV strains, 
respectively, and shared 98.4%-99.4% nt identity and 98.8%-99.4% 
deduced aa identity with the other Chinese strains. The HB-BD 
strain had the highest nt identity (99.4%) and aa identity (99.4%) 
with HKU15-155 and had the lowest nt homology (97.8%) with two 
strains from Thailand, $5011 and S5015L. Further analysis showed 
that the HB-BD strain had one aa substitution (V 43 A), which is 
similar to observations in the HKU15-155 and CHN-AH-2004 
strains. In addition, HB-BD had one further aa mutation (A 308 V). 
To analyse the phylogenetic relationships between the HB-BD 
isolate and the reference strains, a phylogenetic tree was con- 
structed. The representative tree is shown in Figure 3. In the context 


of the complete genome phylogenetic tree, all 29 PDCoV strains 
were separated into the genus Deltacoronavirus, a cluster that was 
distinct from the PEDV and TGEV strains. 

Further analysis demonstrated that the PDCoV isolate, HB-BD, 
was more closely related to other Chinese PDCoV isolates than to 
those isolated from the United States, South Korea and Thailand 
(Figure 3). 


4 | DISCUSSION 


PDCoV is novel swine enteropathogenic coronaviruses that have 
emerged in several pig breeding countries in recent years (Homwong 
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FIGURE 3 _ Phylogenetic tree of the nucleotide sequences of PDCoV isolates based on the complete genome. All the reference sequences 
used in this study were obtained from the GenBank database. The names of the strains, places of isolation and GenBank accession numbers 
are shown in Table 1. The Tree was constructed by the neighbour-joining method with 1,000 bootstrap replications using MEGA 7.0.14 
software. The isolate identified in this study is indicated by black triangle, and the Chinese reference strains are indicated by black dots 


et al., 2016; Madapong et al., 2016; Zhai et al., 2016; Zhang, 2016). 
In China, as PDCoV HKU15-44 and HKU15-155 strains were first 
identified in Hong Kong in 2012 (Woo et al., 2012), PDCoV has 
been reported in some local provinces and the prevalence of PDCoV 
has been investigated (Chen, Zhu, et al., 2015; Dong et al., 2015; 
Song et al., 2015; Wang et al., 2015). In one study, 64 faecal and 
intestinal samples collected from six pig farms in the Shanxi, Guang- 
dong and Hubei provinces in China, were examined for enteropatho- 
gens using RT-PCR (Chen, Zhu, et al., 2015). The prevalence of 
PDCoV infection was 23.4% (15/64). In one survey performed by 
Song et al. (2015), 356 porcine faecal/intestinal samples that were 
collected between November 2012 and March 2015 from diarrhoeic 
pigs on 51 farms in Jiangxi province, China, were examined using 
nested RT-PCR. The result showed that monoinfection with PDCoV 
was 33.71% (120/356) and coinfection with PDCoV and PEDV was 


19.66% (70/356), thus 58.33% (70/120) of the PDCoV-positive sam- 
ples were co-infected with PEDV. In the study performed by Dong 
et al. (2015), 215 intestinal or faecal samples were collected from 
piglets with clinical diarrhoea at various times during 2004-2014 in 
Anhui, Guangxi, Hubei and Jiangsu provinces, of mainland China, for 
the detection of enteropathogens. Among these samples, 14 (6.51%) 
were positive for PDCoV and 50% (7/14) of the 14 PDCoV-positive 
samples were also positive for PEDV. These results showed that 
PDCoV coinfections with other pathogens, especially PEDV, are 
more common in these provinces. However, little is known about 
the prevalence of PDCoV in Hebei province, a major pig-rearing pro- 
vince in China. In the present study, 49 faecal and intestinal samples 
from piglets with diarrhoea on different farms in Hebei province 
were collected, from which PDCoV was detected using RT-PCR. 
Other general enteropathogens, such as PEDV, TGEV and PRoV, 
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were also examined. Five (10.2%) of the 49 samples were positive 
for PDCoV, 31 (63.3%) were positive for PEDV, and two (4.1%) were 
positive for PRoV. Although the prevalence of PDCoV infection was 
lower than that of PEDV, the result indicates that PDCOV infection 
exists in the Hebei pig population. Interestingly, no coinfection with 
PDCoV and PEDV was detected. In addition, all five PDCoV-positive 
samples were negative for all other enteropathogens detected in the 
present study. These findings are somewhat different from observa- 
tions of infection with PDCoV in other provinces and other coun- 
tries. If a greater number of samples were collected, some 
coinfection might have been detected. 

Although the prevalence of PDCoV in some provinces of China 
has been confirmed and the genetic characteristics of some Chinese 
PDCoV strains and those from other countries have been analysed, 
there is little published information about Chinese PDCoV isolates. 
In the present study, using ST cells, the PDCoV HB-BD strain was 
successfully isolated from the intestinal contents of a piglet with 
diarrhoea. After several passages, propagation of PDCoV HB-BD 
was confirmed by the detection of a typical CPE and infectious virus 
titration. The M, N and S gene sequences of HB-BD were also deter- 
mined. To confirm whether the replication of any other swine 
enteric viruses occurred on ST cells, the cell culture-passaged sam- 
ples were detected by RT-PCR and tested negative for PDCoV, 
PEDV, TGEV and PRoV. These results all indicate that the PDCoV 
HB-BD isolate was highly adapted to ST cells. The findings of Hu, 
Jung, et al. (2015) and Dong et al. (2016) indicate that the success 
rate of PDCoV isolation from samples of intestinal contents Is rela- 
tively higher than that from faecal samples. However, even if the 
PDCoV virus is isolated from intestinal contents, the success rate is 
also influenced by other factors, such as virus viability, infectious 
virus titres and other substances in the intestinal contents (Hu, Jung, 
et al., 2015; Dong et al., 2016). In the present study, we attempted 
to isolate the virus from two intestinal content samples and three 
faecal samples that were positive for PDCoV. However, in the ST 
cells, only the PDCoV HB-BD strain was successfully isolated from 
intestinal contents. The success rate of virus isolation is relatively 
low. These results are consistent with those of previous reports 
(Dong et al., 2016; Hu, Jung, et al., 2015). 

To understand the genetic characteristics of the HB-BD isolate, 
the complete genome of the isolate was sequenced, and phyloge- 
netic analysis was performed with known sequences selected from 
the GenBank database. The results showed that all PDCoV strains 
were distinct from PEDV and TGEV, and clustered within the genus 
Deltacoronavirus. The HB-BD strain had a closer relationship with 
other Chinese strains than those from other countries and was 
grouped within the Chinese PDCoV cluster in the phylogenetic trees. 
When compared with the Chinese strains, HB-BD strain showed 
most closely related to the NH strain (emerged in Heilongjiang pro- 
vince) and the HKU15-155 strain (isolated in Hong Kong). Based on 
the S gene, HB-BD shared the highest aa identity (99.2%) with the 
NH strain. In addition, both strains have two aa mutations (L45H 
and Y123H). However, based on the N gene, the HB-BD strain had 
the highest aa identity (99.4%) with HKU15-155 and shared a similar 
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aa substitution (V43A). Further analysis showed that HB-BD had 
four further aa mutations (H149Y, R888T, A894T and G910C) in the 
S gene, one further aa mutation (ASO8V) in the N gene and a unique 
nucleotide mutation (C63T) in the M gene. Studies have shown that 
the S protein of coronaviruses is associated with receptor binding 
and host adaptation (Graham & Baric, 2010; Sato et al., 2011) and 
the N protein participates in viral transcription and replication (Brian 
& Baric, 2005). Further research in HB-BD should be conducted to 
understand whether these mutations can lead to changes in viru- 
lence, pathogenicity and antigenicity of PDCoV. 

In conclusion, the presence of PDCoV in Hebei province was 
confirmed, and the PDCoV strain, HB-BD, was successfully isolated 
and serially passaged in ST cells. The complete genome of HB-BD 
was sequenced and characterized. These results will be valuable for 
further research of PDCoV genetic evolution. The immunogenicity 
and pathogenicity of the HB-BD strain require further study for the 
development of effective diagnostic reagents, assays and potential 


vaccines against newly emerged PDCoV strains. 
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